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ABSTRACT 

Many  important  Fredholm  integral  equations  have  separable  kernels, 
which  are  finite  rank  modifications  of  Volterra  kernels.  This  class  of  kernels 
includes  Green's  functions  for  Sturm-Liouvi 1 le  and  other  two-point  boundary 
value  problems  for  linear  ordinary  differential  operators.  It  is  shown  how  to 
construct  the  Fredholm  determinant,  resolvent  kernel,  and  eigenfunctions  of 
kernels  of  this  class  by  solving  related  Volterra  integral  equations  and  finite 
linear  algebraic  systems.  Applications  to  boundary  value  problems  are  discussed 
and  explicit  formulas  are  given  for  a simple  example.  Analytic  and  numerical 
approximation  procedures  for  more  general  problems  are  indicated. 
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RESOLVENT  KERNELS  OF  GREEN’S  FUNCTION  KERNELS  AND  OTHER 
FINITE  RANK  MODIFICATIONS  OF  FREDHOLM  AND  VOLTERRA  KERNELS 

L.  B.  Rail 

1.  A classical  problem  in  the  theory  of  integral  equations.  The  equation 

1 

(1-1)  y(x)  - a / K(x,t)y(t)dt  = f(x)  , 0 < x < 1 , 

0 = = 

for  the  unknown  function  y(x)  is  called  a linear  integral  equation  of 
second  kind.  Equations  of  this  form  arise  in  the  solution  of  initial  and 
boundary  value  problems  for  ordinary  differential  equations  and  in  other  areas 
of  applied  analysis.  In  (1.1),  the  function  f(x),  the  kernel  K(x,t),  and 
the  parameter  a are  assumed  to  be  given.  In  case  K(x,t)  = 0 for  t > x, 
equation  (1.1)  is  said  to  be  of  Vol terra  type,  and  the  interval  of  integration 
is  actually  0 < t < x ; otherwise,  (1.1)  is  called  an  integral  equation  of 
Fredholm  type. 

A central  problem  in  the  classical  theory  of  linear  integral  equations 
of  second  kind  is  to  determine  the  values  of  A for  which  the  solution  y(x) 
of  (1.1)  exists  and  is  unique,  and  to  express  this  solution  in  the  form 

1 

(1.2)  y(x)  = f(x)  + a / R(x,t;A)f(t)dt  , 0 < x < 1 , 

0 

where  R(x,t;A)  is  called  the  resolvent  kernel  of  K(x,t)  [12].  The  investi- 
gation of  the  unique  solvability  of  equation  (1.1)  can  thus  be  reduced  to  the 
problem  of  existence  and  construction  of  R(x,t;A)  . 

Using  operator  notation,  the  kernel  K(x,t)  may  be  taken  to  define  the 


Sponsored  by  the  United  States  Army  under  Contract  No.  DAAG29-75-C-0024. 


linear  integral  operator  K on  the  space  of  functions  considered.  If  I de- 
notes the  identity  operator,  then  (1.1)  can  be  written  in  the  form 
(1.3)  (I  - xK)y  = f . 

The  solution  y of  (1.3)  exists  and  is  unique  if  the  operator  I - xK  is 
invertible  for  the  given  value  of  x . To  obtain  the  expression  corresponding 
to  (1.2),  the  inverse  of  I - xK  is  represented  in  the  form 


(1.4)  (I  - xK)"1  = I + xR(x)  , 

where  the  resolvent  operator  R(x)  of  K is  the  linear  integral  operator  with 
kernel  R(x,t;x)  . Equation  (1.4)  leads  directly  to  the  relationships 

R(x)  « K + xKR(x)  , 

R(x)  = K + xR(x)K  , 
by  the  definition  of  the  inverse  operator.  These  are  the  so-called  resolvent 
equations.  In  terms  of  the  corresponding  kernels,  equations  (1.5)  become 


(1-5) 


(1.6) 


1 

R(x,t;x)  = K(x,t)  + xf  K(x,s)R(s,t;x)ds  , 

0 

1 

R(x,t;x)  = K(x,t)  + xf  R(x,s;x)K(s,t)ds  . 

0 


(If  K(x,t)  is  a Volterra  kernel,  then  the  intervals  of  integration  in  (1.6) 
reduce  to  t < s < x . ) 


In  the  classical  setting  of  the  theory  of  integral  equations,  one  is 
concerned  with  kernels  which  are  bounded  and  at  least  square-integrable.  The 
representation  (1.4)  allows  one  to  eliminate  the  identity  operator  and  obtain 
relationships  (1.5)  between  linear  integral  operators  with  kernels  (1.6)  of 
this  type.  This  simplifies  the  analysis  considerably,  as  the  identity  operator 
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cannot  be  represented  as  a linear  integral  operator  with  a bounded  kernel 
on  the  spaces  of  continuous  or  square-integrable  functions  [6,  8,  16]  . 


2.  Vol terra  resolvent  kernels.  In  case  that  K(x,t)  is  a Volterra  kernel, 

(1.1)  takes  the  form 

x 

(2-1)  y(x)  - xf  K(x,t)y(t)dt  = f(x)  , 0 < x < 1 . 

0 

It  is  well-known  [0,  11,  12,  16]  that  the  resolvent  kernel  R(x,t;A)  of 
K(x,t)  always  exists  under  the  assumptions  of  the  classical  theory,  and  is 
given  by  the  Neumann  series 

(2.2)  R(x,t; a)  = l An‘Vn)(x,t)  , 0 < t < x < 1 , 

n=l 


where 


(2.3) 


K(1)(x,t)  = K(x,t)  , 

K(n+1)(x,t)  = JX  K(x,s)K(n)(s,t)ds  = /X  K{n)(x,s)K(s,t)ds  , 
t t 


n - 1,  2, 


T 


The  convergence  of  (2.2)  for  all  A with  finite  modulus  is  easy  to  establish 
by  mathematical  induction  for  bounded  kernels  K(x,t)  . If  |K(x,t)|  < M, 

0 < t < x < 1 , then 


|K(n)(x.t) | < M 


(n-1 ) ! 


from  which  the  desired  result  follows.  The  kernels  K^(x,t),  K^(x,t),... 


are  sometimes  called  the  iterated  kernels  of  K(x,t)  [7]  . 

It  will  be  useful  later  to  consider  also  the  transposed  Volterra  Integral 
equation  corresponding  to  (2.1), 

1 

(2.5)  z(t)  - xf  z(x)K(x,t)dx  - g(t)  , 0 $ t i 1 . 

t 
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' 


* 


The  solution  z(t)  of  (2.5)  is  given  in  terms  of  g(t)  and  the  resolvent 
kernel  R(x,t;x)  of  K(x,t)  as 

1 

(2.6)  z ( t ) = g(t)  + a/  g(x)R(x,t;A)dx  , 0 < t < 1 . 

t 


3.  Fredholm  resolvent  kernels.  For  the  more  general  case  of  a Fredholm 
kernel  K(x,t),  an  expression  for  the  resolvent  kernel  is  sought  in  the  form 


(3.1) 


R(x,t;x)  = 


with  the  numerator  and  denominator  having  series  expansions 


(3.2) 


N(x,t;x)  * l An''Kn(x,t) 
n=l 


(3.3) 


a(a)  = 1 + l Cnxn  , 
n=  1 


respectively,  which  converge  for  all  A with  finite  modulus.  The  resolvent 
kernel  R(x,t;x)  will  then  exist  for  all  values  of  X for  which  the 
Fredholm  determinant  a(a)  of  K(x,t)  does  not  vanish.  This  is  analogous 
to  Cramer's  rule  for  the  inversion  of  a finite-dimensional  matrix. 

Formulas  for  the  so-called  associated  kernels  K^(x,t),  K2(x,t),  ... 
of  K(x,t)  [7]  appearing  in  (3.2)  may  be  obtained  by  substituting  (3.1)  - 
(3.3)  into  the  resolvent  equations  (1.6).  This  gives 

^(x.t)  = K(x,t) 


(3.4) 


Kn+1(x,t)  = cnK(x,t)  + / K(x,s)Kn(s,t)ds  « 


= c K(x,t)  + / K (x,s)K(s,t)ds,  n = 1,2,..., 
n o 

a result  which  satisfies  (1.6)  formally,  independently  of  the  values  assigned 
to  c,,  c2,  ...  • For  example,  if  c^  = c2  = •••  = 0,  then  (3.2)  becomes 


4m 


the  Neumann  series  (2.2),  which  does  not  converge  in  general  for  |x| 
large.  In  order  for  (3.2)  and  (3.3)  to  be  entire  functions  of  x , one 
chooses 


(3.5)  cn  = ' n tr  Kn  = ' n ^ Kn(x»x)dx  * n = . 

where  the  quantity  tr  Kn  is  called  the  trace  of  the  kernel  Kn(x,t)  . The 

construction  (3.4)  - (3.5)  of  the  resolvent  kernel  is  in  the  form  given  by 
Lalesco  [11].  Fredholm  [5]  originally  obtained  the  formulas 


tn-l  1 1 


K(s,t)  K(x,s1)  •••  K(x,snl) 


(3.6)  Kn(»,t)  - K(s, ,t)  Msj.s,)  ds, 


‘'Vi'11  K(sn-rsl)'"K(sn-l-sn-l) 


and  the  similar  expressions  corresponding  to  (3.5)  for  cn,  n = 1,2,...  . 

The  satisfaction  of  (3.4)  is  easily  verified  by  mathematical  induction. 

Using  Fredholm's  formulas,  the  convergence  of  the  series  (3.2)  and 
(3.3)  can  be  established  on  the  basis  of  Hadamard's  inequality  for  determinants 
(see  [16]  for  an  elegant  proof)  and  the  ratio  test.  If  |K(x,t)|  $ M,  then 
Hadamard's  inequality  applied  to  (3.6)  yields  the  estimates 


(3.7) 


Mn-12 

| K (x,t)  | < M • L-JJ-  , n = 1,2 

" (n-1)! 


and  similar  bounds  for  Ic^,  | c^  | . . . . . The  rate  of  convergence  which  can 

be  predicted  for  the  series  (3.2)  and  (3.3)  on  the  basis  of  (3.7)  is,  of 
course,  much  slower  than  that  given  by  the  estimates  (2.4)  applied  to  the 


Neumann  series  expansion  (2.2)  of  a Volterra  resolvent  kernel.  For  example, 
if  given  Volterra  and  Fredholm  kernels  are  bounded  in  absolute  value  by  M , 
then  for  n * 100,  the  right-hand  side  of  (3.7)  exceeds  that  of  (2.4)  by 
a factor  of  10100  . From  a computational  point  of  view,  the  relationships 
(3.4)  - (3.5)  would  appear  to  be  preferable  to  the  equivalent  expressions 
(3.6)  involving  determinants.  It  has  been  observed,  however,  that  the  for- 
mulas corresponding  to  (3.4)  - (3.5)  for  finite-dimensional  matrices  are 
unstable  numerically  [4]  . 

Other  important  relationships  which  follow  from  the  formulation  (3.1) 
of  the  resolvent  kernel  and  the  resolvent  equations  (1.6)  are 


(3.8) 


1 

N(x,t; a)  = a( a)K(x, t)  + A / K(x,s)N(s,t;A)ds  , 


1 

N(x,t;A)  = a( A)K(x,t)  + A / N(x,s; A)K(s,t)ds  . 

0 


These  are  immediately  evident  for  a(a)  + 0 . In  the  framework  of  the 
classical  theory,  they  can  also  be  extended  to  the  case  that  A = A*  is 
an  eigenvalue  of  the  kernel  K(x,t);  that  is,  a( a*)  = 0 . Assuming  that 
N(x,t;A*)  does  not  vanish  identically,  a point  U.t)  in  the  square 
0 < x,t  < 1 exists  such  that  the  functions 


(3.9) 


y*  (x)  = N(x,t;a*),  0 < x < 1 

z * (t)  = N(s,t; A*),  0 < t < 1 


\ 

are  also  not  identically  zero.  Furthermore,  y*(x) 
geneous  integral  equation 


» 


satisfies  the  homo 


8- 


I 


(3.10)  y*(x)  = A*/  K(x,t)y*(t)dt  , 0 < x < 1 , 

0 

and  is  said  to  be  a (right)  eigenfunction  of  K(x,t)  corresponding  to 
the  eigenvalue  A*  . Similarly,  the  function  z*(t)  satisfies  the 

transposed  homogeneous  integral  equation 

(3.11)  z*(t)  = A*/  z*(x)K(x,t)dx  , 0 < t < 1 , 

0 

and  is  called  a left  eigenfunction  of  K(x,t)  corresponding  to  A*  . 


4.  Symmetric  separable  kernels.  Attention  will  now  be  devoted  to  the 
construction  of  resolvent  kernels  of  a special  class  of  Fredholm  kernels. 
In  the  symmetric  case,  a kernel  of  the  form 


(4.1) 


G(x,t)  - 


u(t)v(x) , 0 < t * x s 1 , 


u(x)v(t) , 0 < x < t < 1 


will  be  called  a simple  separable  kernel.  It  is  assumed  that  the  functions 
u(x)  and  v(x)  are  linearly  independent;  otherwise,  G(x,t)  would  be  a 
degenerate  kernel  of  rank  one  [8,  pp.  37-40].  In  general,  a symmetric 
separable  kernel  is  a finite  sum  of  linearly  independent  kernels  (4.1)  . 

Before  dealing  with  the  general  case,  the  resolvent  kernel  of  the 
symmetric  simple  separable  kernel  (4.1)  will  be  constructed.  To  do  this, 
the  Fredholm  integral  equation 


(4.2) 


y(x)  - xj  G(x,t)y(t)dt  = f(x),  0 < x < 1 , 

0 


will  be  solved.  This  is  essentially  the  approach  used  by  Drukarev  [3], 
Brysk  [2],  and  Aalto  [1].  Using  the  definition  (4.1)  of  G(x,t),  equation 

(4.2)  may  be  written  as 

x 1 

(4.3)  y(x)  - a / u(t)v(x)y(t)dt  = f(x)  + a / u(x)v( t)y(t)dt  . 

0 x 


Adding  the  quantity  a / u(x)v(t)y(t)dt  to  both  sides  of  (4.3)  gives 

0 


(4.4) 


y(x)  - A / K(x,t)y( t ) d t = f(x)  + Acu(x)  , 
0 


where 


.. 


1 

(4.5)  c = / v(t)y(t)dt 

0 

is  to  be  determined,  and  K(x,t)  is  the  Volterra  kernel 


(4.6)  K(x,t)  = u(t)v(x)  - u(x)v(t),  0 < t < x < 1 . 


The  system  of  equations  (4.4)  - (4.5)  is  easily  seen  to  be  equivalent  to 
the  original  integral  equation  (4.2). 

As  was  shown  in  §2,  the  kernel  K(x,t)  has  the  Volterra  resolvent 
kernel  R(x,t;A)  given  by  (2.2)  for  all  A . Define 


r x 

F(x)  = f ( x)  + A / R(x,t;x)f(t)dt  , 

0 

(4'7) 

U(x)  = u(x)  + A / R(x,t;A)u(t)dt  , 

^ 0 

where  the  dependence  of  F(x)  and  U(x)  on  A has  been  suppressed  for 
clarity  of  notation.  From  (4.4), 


(4.8)  y(x)  = F(x)  + AcU(x)  , 


and  from  (4.5) 

1 1 

(4.9)  c = / v(t)F(t)dt  + Ac  / v(t)U(t)dt  . 

0 0 

Thus,  (4.9)  has  a unique  solution  for  c if 

1 

(4.10)  A(A)  = 1 - a/  v(t)U(t)dt  ^ 0 , 

0 


in  which  case 
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(4.11)  c = } v(t)F(t)dt  . 

a(  X)  0 

The  expression  (4.11)  for  c may  be  written  in  terms  of  f(t)  by  intro- 


ducing the  function 


(4.12) 


V(t)  = v(t)  + A / v(x)R(x,t;A)dx  . 
t 


It  then  follows  from  (4.7)  and  (4.11)  that 


(4.13) 


and,  from  (4.8), 


c « J V(t)f(t)dt  , 
M a)  0 


(4.14) 


y(x)  = f(x)  + a/  R(x,t;A)f(t)dt  + / V(t)f(t)dt  . 

0 4lA;  0 


By  definition,  the  resolvent  kernel  r(x,t;x)  of  G(x,t)  can  be  obtained 


directly  from  (4.14)  as 


(4.15)  r(x,t;x)  = 


R(x,t;x)  + -TT  U(x)V(t),  0 < t £ x =<  1 , 


U(x)V(t), 


0 < x 5 t 5 1 , 


provided,  of  course,  that  A( X ) f 0 . 

Another  expression  for  r(x,t;A)  can  be  obtained  by  making  use  of 
the  fact  that  the  resolvent  kernel  of  a symmetric  kernel  G(x,t)  must  also 
be  symmetric  in  x and  t;  that  is,  r(x,t;A)  = r(t,x;x).  This  gives 


U(t)V(x) , 


0 £ t < X < 1 , 


(4.16)  r(x,t;x)  * { AVAJ  i , t „ . 

| R(t,x;x)  + U(t)V(x),  0<x<t<l  , 

which,  when  compared  with  (4.15),  yields 


(4.17)  a( x)R(x,t; X)  = U( t)V(x)  - U(x)V(t),  0 < t < x < 1 , 

and  finally. 


r 


0^* 


0 < t < x < 1 , 

0 < x < t < 1 . 

These  results  give  rise  to  the  following  observations. 

Remark  4.1.  Equation  (4.18)  shows  that  if  the  resolvent  kernel  of 
a symmetric  simple  separable  kernel  (4.1)  exists,  then  it  is  also  a symmetric 
simple  separable  kernel. 

Remark  4.2.  The  expression  (4.18)  for  the  resolvent  kernel 
r(x,t;x)  of  G(x,t)  can  be  obtained  without  finding  the  Volterra  resolvent 
kernel  R(x,t;x)  of  K(x,t)  explicitly;  one  need  only  solve  the  Volterra 
integral  equation  (2.1)  with  f(x)  = u(x)  for  y(x)  = U(x),  and  the 
transposed  Volterra  integral  equation  (2.5)  with  g(t)  = v(t)  for 
z(t)  = V(t)  . 

Remark  4.3.  An  equivalent  expression  for  the  Fredholm  determinant 
A(x)  of  G(x,t)  is 

1 

(4.19)  A ( X ) = 1 - x/  V(x)u(x)dx  , 

0 

which  can  be  obtained  directly  from  (4.10)  by  interchange  of  the  order  of 
integration  and  use  of  the  definition  (4.12)  of  V(t).  It  is  easy  to  show 
that  the  above  expressions  give  expansions  of  a(x)  and  N(x,t;x)  in 
powers  of  x with  coefficients  satisfying  the  relationships  (3.4)  - (3.5) 

[1,  13].  The  advantage  of  the  present  approach  is  that  the  rate  of 
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(4.18) 


r(x,t;x)  = 


1 

a(x) 

1 

a(x) 


U( t)V(x) , 
U(x)V(t) , 


convergence  of  the  series  for  a(x),  U(x)  = U(x;x),  and  V(t)  = V(t;x) 
can  be  predicted  on  the  basis  of  (2.4)  rather  than  (3.7). 

Remark  4.4.  The  definitions  of  the  Vol terra  kernel  K(x,t)  and 
its  resolvent  kernel  R(x,t;A)  extend  to  the  entire  square  0 < x,  t < 1, 
with 

(4.20)  K(x,t)  = 0,  R(x,t;A)  =0,  0 < x < t < 1 . 

This  will  result  in  continuous  kernels  for  u(x),  v(t)  continuous,  as 
K(x,x)  = R(x,x;A)  = 0,  0 < x < 1 . In  terms  of  these  extended  kernels, 
one  may  write 

f G(x,t)  = K(x,t)  + u(x)v(t)  , 

(4.21)  { . 

1 r(x,t;X)  = R(x,t;x)  + — U(x)V(t)  . 
v A(x) 

Thus,  the  kernel  G(x,t)  given  by  (4.1)  is  the  sum  of  the  Volterra 
kernel  K(x,t)  and  the  degenerate  kernel  u(x)v(t)  of  rank  one.  If 

a(x)  + 0,  then  its  resolvent  kernel  r(x,t;x)  exists  and  is  also  a 
rank  one  modification  of  the  Volterra  resolvent  kernel  R(x,t;x)  of 
K(x,t).  By  symmetry,  one  also  has 

G(x,t)  = K(t,x)  + u ( t ) v ( x)  , 

(4.22)  / 

r(x,t;x)  * R(t,x;A)  + — U(t)V(x)  . 

v a(x) 

In  some  applications,  the  expressions  (4.21)  or  (4.22)  for  G(x,t)  arise 
more  naturally  than  (4.1). 

The  next  case  to  be  examined  is  that  x = A*  is  an  eigenvalue  of 
G(x,t);  that  is,  a(x*)  * 0 . Noting  that  U(x)  = U(x;x)  and 
V(t)  * V(t;x)  depend  on  x , define 


I 
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U*(x)  = U(x;x*) , 0 < x < 1 , 

(4.23)  ( 

1 V*(t)  = V(t;x*),  0 < t < 1 . 

From  (4.17),  which  can  be  extended  easily  to  x = x*,  it  follows  that 
the  functions  U*(x)  and  V*(x)  are  linearly  dependent.  Hence,  for 
x = x*,  the  Fredholm  numerator  N(x,t;X)  of  r(x,t;x)  becomes 

(4.24)  N(x,t;X*)  = aU*(x)V*(t)  , 

where  a f 0 is  some  constant.  Thus,  by  (3.9)  - (3.11),  U*(x)  will  be 

a right  eigenfunction  of  G(x,t)  corresponding  to  x*,  and  V*(t)  is  a 
corresponding  left  eigenfunction.  As  G(x,t)  is  symmetric,  the  distinction 
between  left  and  right  eigenfunctions  is  inconsequential.  However,  it  will 
be  shown  later  that  a similar  approach  gives  these  as  distinct  functions  in 
the  nonsymnetric  case. 

The  implications  of  the  above  results  for  the  solvability  of  the 
integral  equation  (4.2)  may  be  summed  up  in  the  following  familiar  language. 

Theorem  4.1  (Fredholm  Alternative).  If  the  transposed  homogeneous 
integral  equation 

1 

(4.25)  z ( t)  - X/  z(x)G(x,t)dx  =0,  0 5 t < 1 , 

0 

has  only  the  trivial  solution  z(t)  = 0,  0 < t < 1 , then  equation  (4.2) 
has  the  unique  solution 

(4.26)  y(x)  = F(x)  + / v(t)F(t)dt,  0 < x < 1 . 

a(x)  o 

On  the  other  hand,  if  (4.25)  has  the  nontrivial  solution  z(t)  = V*(t), 
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then  (4.2)  has  a solution  only  if  f(x)  is  orthogonal  to  V*(x),  that 


(4.27) 


1 

/ V*(t)f(t)dt  = 0 
0 


If  (4.27)  is  satisfied,  then  (4.2)  has  the  solutions 


(4.28) 


y(x)  = F(x)  + aU*(x)  , 


where  a is  arbitrary,  and  y(x)  = U*(x)  is  a nontrivial  solution  of 
the  homogeneous  integral  equation 


(4.29)  y(x)  - \J  G(x,t)y(t)dt  =0,  0 5 x < 1 . 

0 

As  derived  above,  of  course.  Theorem  4.1  applies  only  to  simple 
kernels  (4.1).  The  same  technique,  however,  applies  to  general  symmetric 
separable  kernels.  Suppose  that 


(4.30)  G(x,t)  = l G.(x.t)  «< 


l U.(t)v.(x),  0 < t 5 x < 1 , 

j=l  J J 


l U.(x)v.(t),  0 < x 5 t < 1 

j-1  J J 


Then,  the  integral  equation  (4.2)  with  the  syirmetric  separable  kernel 

(4.30)  may  be  reduced  to 

(4.31)  y(x)  - xj  K(x,t)y(t)dt  = f(x)  + X l c.  u.(x)  , 

0 j=l  J J 

where  the  numbers 


' 


(4.32) 


ci  = / v.(t)y(t)dt,  j = 1,2 n , 

J n J 


are  to  be  determined,  and  K(x,t)  is  the  Vol terra  kernel 

(4.33)  K(x,t)  = l [u.(t)v.(x)  - U.(x)v.(t)],  0 < t < x < 1 . 

j=l  J J J J 

As  before,  let  R(x,t;X)  denote  the  Volterra  resolvent  kernel  of  K(x,t), 


and  defjpe 


(4.34) 


F(x)  = f(x)  + x/  R(x,t;X)f(t)dt  , 
0 


U,(x)  = u,(x)  + x/  R(x,t;X)u.(t)dt,  i = 1,2,. ...n 
1 0 1 


Equation  (4.31)  is  then  equivalent  to 


(4.35) 


y(x)  = F(x)  + X l U. (x)  . 

i=l  1 1 


Multiplying  (4.35)  by  v^(x),  v^(x),  ....  v^(x)  in  turn  and  integrating 

with  respect  to  x from  0 to  1 gives  the  system  of  equations 

n 1 1 

(4.36)  c.  - X l c J v.(x)U,(x)dx  = / v.(x)F(x)dx,  i = 1,2,... ,n  , 

1 j=l  J 0 1 J 0 1 

for  the  unknowns  c^,  ...  , cp  defined  by  (4.32). 

Remark  4.5.  It  follows  from  (4.32),  (4.35),  and  (4.36)  that  the 
solution  of  the  integral  equation  (4.2)  with  the  symmetric  separable  kernel 
(4.30)  is  equivalent  to  solving  the  Volterra  integral  equation  (2.1)  with 
kernel  (4.33)  and  right-hand  sides  f(x),  u^x),...,  un(x)  for 


2 

F(x),  U,(x),...,  Un ( x ) , forming  the  n + n coefficients  of  the  system 

(4.36) ,  and  then  solving  this  system  of  n linear  equations  for  the  n 
unknowns  c^,  c2,...,Cn  . 

It  is  also  possible  to  obtain  explicit  expressions  for  the  Fredholm 
determinant  a(a)  of  G(x,t)  and  its  resolvent  kernel  r(x,t;A)  in  terms 
of  the  determinant  and  the  inverse  of  the  coefficient  matrix  of  the  system 

(4.36) .  It  will  be  convenient  to  introduce  the  functions 

1 

(4.37)  V - ( t ) = v . ( t)  + xj  v.(x)R(x,t;x)dx,  i = 1,2,.. .,n  , 

1 1 t 

which  are  the  solutions  of  the  transposed  Volterra  integral  equations  (2.5) 
with  right-hand  sides  v^(t),  v^tt),  ...  , vn ( t ) . Let  the  coefficients 

1 1 

(4.38)  = / v.(x)U.(x)dx  = / V.(t)u,(t)dt,  i , j = 1,2,. ...n  , 

ij  0 J 0 J 

define  the  n x n matrix  A = (a-)  • Then,  the  coefficient  matrix  of 

' J 

the  system  (4.36)  has  the  form 

(4.39)  I - xA  - (6ij  - Aoij)  , 

where  6^  is  the  Kronecker  delta  = 0 if  i f j,  6^  = 1),  and 

I = ’s  the  n x n identity  matrix.  If  the  determinant 

(4.40)  &( a)  = det  (I  - aA) 

does  not  vanish,  then  the  inverse  of  the  matrix  I - aA  exists,  and 
may  be  written  as 


(4.41) 


(I  - AA)' 


a(  x) 


B(A) 


by  the  use  of  Cramer's  rule.  In  this  case,  the  system  (4.36)  has  the 
unique  solutions 


(4.42) 


1 n 1 

c.  = -7-r  l 6i i ( X)  / v,(x)F(x)dx  = 
1 AU)  j=i  ^ o 3 


i(IT  ,1,  6u(1)  l0  Vt)f(t)dt 


for  i = l,2,...,n  . Thus,  the  integral  equation  (4.2)  will  have  the 
unique  solution  y(x)  given  by  (4.35),  which  may  be  written  in  the  form 


x 

y(x)  = f(x)  + \j  R(x,t;x)f(t)dt  + 

0 

(4.43) 

n n B,,(x)  1 

♦ » l l U,(x)  ,J  / V,(t)f(t)dt  . 

“mTT  0 J 


From  (4.43),  it  is  possible  to  derive  several  expressions  for  the 
resolvent  kernel  r(x,t;*)  of  the  syrnnetric  separable  kernel  G(x,t), 
provided  that  the  Fredholm  determinant  a(x)  of  G(x,t)  is  nonzero. 
Extending  the  Volterra  kernel  K(x,t)  and  its  resolvent  kernel  R(x,t;x) 
to  the  entire  square  0 < x,  t < 1 as  before,  one  may  write 


(4.44) 


G(x,t)  = K(x,t)  + l u.(x)v,(t)  , 

j=l  J J 

r(x,t;x)  = R(x,t;x)  + — r l £ U . (x) B.  .(x)V .(t)  . 

A \ A / •)  = ] j=]  ' 'J  J 


Remark  4.6.  The  symmetric  separable  kernel  (4.30)  is  a rank  n 
modification  of  the  Volterra  kernel  (4.33).  If  A(x)  f 0,  the  resolvent 
kernel  r(x,t;x)  is  likewise  a rank  n modification  of  the  Volterra 
resolvent  kernel  R(x,t;x)  of  K(x,t)  . 

Integral  operators  corresponding  to  degenerate  kernels  of  rank  n 
are  sometimes  called  n-term  dyads  [6].  By  symmetry. 


(4.45) 


G(x,t)  = K(t,x)  + l u.(t)v.(x)  , 

j=l  J J 

. n n 

r(x,t;X)  = R(t,x;X)  + l £ lUt)B.  ,(X)V.(x)  . 

a(a)  i=!  j=l  1 J 


Comparison  of  (4.45)  and  (4.44)  gives 


n n 


(4.46)  A(A)R(x,t;A)  = £ £ B„(x)[lUt)V,(x)  - MxjV^t)] 

i=l  j=l  1 J 1 J 

and  thus  r(x,t;x)  may  also  be  written  in  the  form 


(4.47)  r(x,t;x)  * 


n n Bn(x) 

[ l -4“r~  U.(t)V.(x),  0 < t < x < 1 , 

i*l  a(x)  1 J 

n n 8n(x) 

l l -^Vx^t),  0 ^ ^ t ^ 1 • 

1=1  j«l  &U)  1 J 


mm 


Remark  4.7.  If  it  exists,  the  Fredholm  resolvent  kernel 
r(x,t;x)  of  the  symmetric  separable  kernel  (4.30)  is  likewise  the  sum  of 
n symmetric  simple  separable  kernels. 

For  n > 1,  there  are  a number  of  ways  in  which  the  kernel  (4.47) 
may  be  written  in  the  form  (4.30).  For  example,  defining  the  functions 
r n 


(4.48)  { 


M*)  = l Ui(x)B..(x),  j = 1 
J i=1  i ij 


9 2 «...  9 n 9 


^(t)  = l i(A)V.(  t), 
1 j-1  1J  J 


i “ 1 ,2,. . . ,n  , 


one  obtains  the  equivalent  representations 

n 

I 

j-1 


' aTxT  I 


(4.49)  r(x,t;x)  = \ 


1 


V. 


HIT  J,  rj(x)vj(t) 


0 < t < x < 1 , 


0 < x < t < 1 , 


and 


^ i " 


HIT  <L 


(4.50)  r(x,t;x)  = l 


i = l 
n 


V. 


HIT  j,  uiW‘i(t> 


0 < t < x < 1 , 


0 < x < t < 1 , 


for  the  resolvent  kernel  of  G(x,t). 

it 

Suppose  now  that  x = x is  an  eigenvalue  of  the  kernel  (4.30), 

★ 

that  is,  a(x  ) = 0.  By  the  same  method  as  used  in  the  nonhomogeneous  case, 
the  homogeneous  equation  (4.29)  can  be  shown  to  be  equivalent  to  the  system 


r y(x)  = a*  l c .u . ( x) 
i = 1 


(4.51) 


V. 


c.“A  £ oijjC.-O,  i - 1 , 2 , . . . , n , 

j = l J J 


corresponding  to  F(x)  = 0 in  (4.35)  and  (4.36).  As  A ( A* ) = 0,  the 
homogeneous  system 

(4.52)  (I  - a*A)c  = 0 

has  m < n linearly  independent  solutions 


(4.53)  c(k)  * (cik),cJk),...,clk))T, 


k l)2j...  jin  . 


'1  ,v7  ’-"'-n 

Corresponding  to  these  solutions,  which  are  right  eigenvectors  of  the  matrix 

A = (a..)  corresponding  to  the  reciprocal  eigenvalue  A , one  obtains  m 
* u 

linearly  independent  eigenfunctions 
(4.54) 


y£(x)  = A * l c|k)U.(x) 
K i=l  1 ’ 


k = 1 ,2,. . . ,n 


. ★ 
of  the  kernel  G(x,t)  corresponding  to  the  eigenvalue  a . Similarly, 

the  transposed  homogeneous  equation  (4.25)  for  a = A*  has  m linearly 

independent  solutions 


(4.55) 


zf(t)  - X*  l d(k)V-(t), 
x i=l  J J 


k ~ 1 > 2 f • . • * m t 


corresponding  to  the  m linearly  independent  solutions 


(4.56)  d(k)  = (djk),d£k) d^k)), 

of  the  transposed  homogeneous  system 

(4.57)  d(  I - a*A)  * 0 
that  is, 

* n 


k * 1,2,..* ,m  , 


Jj 


L di° 


1*1 


1j 


0, 


n . 
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S(“/V 


(4.58) 


j * 1,2,..., 


It  follows  from  (4.36)  that  the  nonhomogeneous  integral  equation  (4.2) 
will  have  no  solutions  unless  the  orthogonality  conditions 

n m 1 1 n / . x * 

(4.59)  l d . / v.(x)F(x)dx  = / l d\*\(t)f(t)dt  = / z.(t)f(t)dt  - 0 

j=l  J 0 J 0 j=l  J J 0 

hold  for  k = 1,2 m;  that  is,  the  right-hand  side  of  (4.2)  must  be 

orthogonal  to  all  solutions  of  the  transposed  homogeneous  equation  (4.25) 

with  x = x*.  If  (4.59)  holds,  then  (4.2)  is  satisfied  by  the  family  of 

sol utions 

* x * m * 

(4.60)  y(x)  = f(x)  + x / R(x,t;x  )f(t)dt  + l Yky. (x)  , 

0 k=l  K K 

with  Y]  *^2’* ‘ • ,Ym  drbitrary* 

The  above  results  can  be  stated  as  the  corresponding  generalization 
of  the  Fredholm  Alternative  Theorem  to  kernels  (4.30)  with  n >1.  Another 
method  fcr  computing  resolvent  kernels  of  separable  kernels  (symmetric 
or  not)  will  be  indicated  in  a later  section. 
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5.  Applications  to  boundary  value  problems.  Symmetric  separable  kernels 
appear  frequently  as  Green's  functions  for  two-point  boundary  value  problems 
for  ordinary  differential  operators  [17].  More  precisely,  suppose  that 
L[  ] is  a linear  ordinary  differential  operator,  and  a solution  y(x) 


of  the  differential  equation 
(5.1) 


L[y(x)]  = h ( x ) 


is  sought  which,  together  with  its  derivatives  of  lower  order  than  the 
order  of  L[  ],  satisfies  given  conditions  at  x = 0 and  x - 1.  If 
this  boundary  value  problem  ha*  a unique  solution  y(x)  which  can  be 
represented  as 


(5.2) 


y(x)  = / G(x,t)h(t)dt 
0 


for  all  functions  h(x)  from  some  class  such  as  continuous  functions, 
then  G(x,t)  is  said  to  be  the  Green's  function  for  L[  ] correspond- 
ing to  the  given  boundary  conditions.  If 


(5.3) 


h ( x ) = Ay(x)  + g ( x ) , 


then  the  boundary  value  problem  for  the  differential  equation  (5.1)  is 
equivalent  to  the  Fredholm  integral  equation 


(5.4) 


y(x)  - x j G(x,t)y( t)dt  * f(x),  0 < x < 1 , 

0 


where 


(5.5) 


f(x)  = / G(x,t)g(t)dt,  0 < x < 1 
0 


Thus,  the  techniques  of  §4  apply  to  the  solution  of  (5.4)  if  G(x,t) 
is  a symmetric  separable  kernel,  and  to  finding  eigenvalues  and  eigen- 
functions in  the  homogeneous  case  g(x)  0. 


In  many  cases,  it  is  possible  to  proceed  directly  from  the 
boundary  value  problem  to  representations  for  the  Green's  function,  its 
Fredholm  determinant,  resolvent  kernel,  and  eigenfunctions.  For  example, 
consider  the  Sturm-Liouvi lie  operator  [17] 


(5.6) 


L[y(x) ] = -(p(x)y' (x) ) 


with  p(x)  > 0,  subject  to  the  boundary  conditions 


(5.7) 


ap(0)y'(0l  + by(0)  = 0 , 
cp(l)y'(l)  + dy ( 1 ) = 0 . 

A simple  way  to  find  the  Green's  function  for  this  problem  is 


to  start  directly  from  the  differential  equation 


(5.8) 

One  integration  gives 


(p(x)y’(x))  = -h(x)  . 


(5.9) 


p(x)y’(x)  = p(0)y'(0)  - / h(t)dt 

0 


Dividing  (5.9)  by  p(x)  and  integrating  again  yields 

X X s 

(5.10)  y(x)  * y(0)  + p(0)y’(0)  / ^- - / / dtds 


By  defining 


(5.11) 


F<’>  ‘l  sm  • 


and  noting  that  change  of  order  of  integration  results  in 

X S x t x X 

(5.12)  £ / ^jdtds  = / / ^fy|h(t)dt  = £ [F(x)  - F(t)]h(t)dt  , 


one  obtains 


(5.13)  y(x)  = y(0)  + p(0)y' (O)F(x)  - / [F(x)  - F(t)]h(t)dt 

0 


from  (5.10).  The  boundary  conditions  (5.7)  will  now  be  used  to  express 
y ( 0 ) and  p(0)y'(0)  in  terms  of  integral  transforms  of  h(x).  From 
(5.9)  and  (5.13), 


(5.14) 


p(Uy'O)  = P(0)y  ’ ( 0) 


1 

/ h(t)dt  , 
0 


y(D  a y(o)  + 


l 


p(0)y' (0)F(1 ) - / 
0 


[F(l ) - F(t)]h(t)dt  . 


Multiplying  the  first  equation  of  (5.14)  by  c,  the  second  by  d,  and 
adding  the  results  gives 

1 

(5.15)  0 = dy(0)  + Yp(0)y'(0)  - / [y  - dF(t)]h(t)dt 

0 


by  the  second  of  the  boundary  conditions  (5.7),  where 

(5.16)  y - c + dF( 1 ) . 

Now,  multiplying  (5.15)  by  a and  using  the  first  boundary  condition  of 
(5.7)  results  in 

1 

(5.17)  0 = (ad  - bY)y(O)  - / a[y  - dF(t)]h(t)dt  , 

0 

or,  if 

(5.18)  6 = ad  - by  t 0 , 
then 

1 , 

(5.19)  y(0)  = / | [y  - dF(t)]h(t)dt  . 

0 6 

Similarly,  multiplication  of  (5.15)  by  b and  use  of  the  first  boundary 
condition  gives 


(5.20) 


! 


p( 0)y ' (0)  = - / j [y  - dF(t)]h(t)dt  . 

0 

Equations  (5.19)  and  (5.20)  may  be  substituted  into  (5.13)  to  obtain 
x 1 1 

(5.21)  y(x)  = - / [ F ( x ) - F(t)]h(t)dt  + / j [a  - bF(x)][y  -dF(t)]h(t)dt  . 

0 0 6 

Comparison  of  (5.21)  with  (5.2)  gives  the  following  results. 

Remark  5,1.  If  6^0,  then  the  Green's  function  G(x,t)  of 
the  Sturm-Liouvi 1 le  problem  (5.6)-(5.7)  is  the  rank  one  modification 

(5.22)  G(x,t)  = K(x,t)  + u(x)v(t) 
of  the  Vol terra  kernel 

(5.23)  K(x.t)  = - [ F ( x ) - F(t)],  0 < t < x < 1 , 

with 


u(x)  = j [a  - bF(x)]  , 


(5.24) 


v(t)  = y - dF( t)  . 

Remark  5.2.  If  0,  then  the  Green's  function  (5.22)  is  the 
symmetric  simple  separable  kernel 


(5.25)  G(x,t)  = 


-1  [a  - bF(t)][y  - dF(x)], 
l [a  - bF(x)][y  - dF(t)]. 


0 < t < x < 1 , 

0 < x < t < 1 . 


This  follows  directly  from  (5.22)-(5.24)  and  the  simple  calculation 
(5.26)  [a  - bF(x)][y  - dF(t)]  - 6F(x)  + 6F(t)  = [a  - bF(t)][y  -dF(x)]  . 

Thus,  the  evaluation  of  the  Fredholm  determinant  a(a)  of  G(x,t) 
and  the  calculation  of  its  resolvent  kernel  and  eigenfunctions  depends 
only  on  being  able  to  calculate  the  functions  U(x),V(t)  either  in  terms 
of  the  Volterra  resolvent  kernel  R(x,t;x)  of  the  kernel  (5.23),  or 
directly  by  solving  the  Volterra  integral  equations  cited  in  Remark  4.2. 
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For  certain  simple  examples,  it  is  possible  to  give  explicit 
formulas  for  these  results.  Taking  p(x)  s 1 gives 

(5.27)  L[y(x) ] = -y"(x)  . 
and  hence 

(5.28)  F(x)  = x,  K(x,t)  = -(x  - t)  , 
from  which 

(5.29)  R(x,t;X)  = - — n*^  ~ ^ . 

The  boundary  conditions 

(5.30)  y(0)  = 0,  y ' ( 1 ) = 0 

correspond  to  (5.7)  with  a=0,  b=l,c=l,d=0,  and  hence  y * 1, 

6 = -1 , and  thus 

(5.31)  u(x)  = x,  v(t)  = 1 . 

It  follows  that  the  Green's  function  G(x,t)  for  the  differential  operator 
(5.27)  with  boundary  conditions  (5.30)  is 

(5.32)  G(x.t)  = -(x  - t)+  + x , 


where  (x  - t)+  = 0 

(5.33) 


for  0 < x < t 5 1 , 


G(x,t) 


t. 


x. 


or 

0 < t < x < 1 
0 < x < t < 1 


Using  (5.29)  and  (5.31),  one  obtains 

(5.34)  U(x)  = , V(t)  = cos/x(l  - t)  , 


from  which  the  Fredholm  determinant  A ( A ) of  G(x,t)  is  seen  to  be 

(5.35)  A(X)  * cos/x  . 

If  a(x)  t 0,  then  the  resolvent  kernel  r(x,t;x)  of  G(x,t)  may  be 
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written  as 


(5.36) 


r(x,t;x) 


sin^x  - t)  + 

7T 


. sin/Xx  cos/x(i  . t) 

+ TTlosh 


corresponding  to  (5.32),  or  as  the  symmetric  simple  separable  kernel 

0 < t < x < 1 , 


(5.37) 


r(x,t;x)  = 


' sin/xt  cos^Q  - x) 
vT  cos  A1 


sin/xx  cos/x(l  - t) 
TTcosA 


0 < x < t < 1 , 


of  the  same  form  as  (5.33).  It  also  follows  immediately  from  (5.35) 
that  the  eigenvalues  of  G(x,t)  are 


(5.38) 


'n  ■[(£=-rLH  • " = '.2.3.. 


(5.40) 


and  the  corresponding  eigenfunctions  are  proportional  to 
(5.39)  /Tu*(x)  = sin  , n*  1,2,3 

It  is  also  possible  to  write  down  explicit  formulas  for  the  more 
general  second-order  boundary  value  problem 

L [y ( x ) ] =-y"(x) 
ay'(0)  + by(0)  = 0 , 
cy'(l)  + dy(l ) = 0 , 
by  the  use  of  (5.24),  (5.28),  and  (5.29).  The  Green's  function  G(x,t) 
is  obtained  from 

u(x)  = 1 [a  - bx]  , 
v(t)  = c + d(l  - t)  , 
provided  6 = ad  - b ( c + d)  f 0,  the  resolvent  kernel  r(x,t;x)  of 
G(x,t)  from 


(5.41) 


U(x)  * jj-  [a  cos/xx  - sin/xx]  , 

V(t)  * c cosAO  - t)  + ^ sin/Kl  - t)  , 


(5.42) 


-29- 


and  the  Fredholm  determinant  is 


(5.43)  [(ad  - be) cos ''T  - (ac^T  + ^sin-'T]  . 

The  eigenvalues  of  G(x,t)  may  thus  be  found  by  solving  the  simple 
transcendental  equation 

(5.44) 

for  and  the  results  substituted  into  (5.42)  to  obtain  the 

corresponding  eigenfunctions.  To  simplify  these  expressions,  one  may 
introduce  the  angles 

(5.45) 


tan /a  _ jid  - be 
/T  Aac  + bd 


which  gives 


(5.46)* 


$ = 4> ( A ) = tan  ^ , e = 0(a)  = tan  ^ , 


U(x)  = 


V(t) 


cos(/ax  + 4>)  , 


sin[/T(l  - t)  + e]  , 


and 

(5.47) 

Thus, 


A(X) 


- I cos(/a  ♦ 4 ♦ 0)  . 


(5.48)  r(x,t;A)  = S, 


( cos(/At  -a  <t> ) s i n [ i/T ( 1 - x)  + oj 
/acos(/T  + 41  + 0) 


0 < t < x < 1 
< x < t < 1 


provided,  of  course,  that  b{\)  t 0. 

From  (5.47),  the  eigenvalues  ApAg****  of  must  satisfy 

(5.49)  /rn  + »(An)  + 0(An)  * (2n-  j 1 ) 7T  , n * 1,2 


1 


¥ 


or 

, (ad  - bc)/T  / 

(5.50)  A «■  cot-'  vc , bd  ■ 

which  is  equivalent  to  (5.44). 

The  method  given  above  extends  readily  to  two-point  boundary  value 
problems  of  arbitrary  order.  In  general,  the  Volterra  kernel  K(x,t)  is 
obtained  by  integrating  (5.1)  as  an  initial  value  problem  with  zero 
initial  conditions;  for  example,  for 

(5.51)  L[y( x) ] = -yiv(x)  , 
one  obtains 

(5.52)  K(x,t)  = - , 0 < t < x < 1 . 

The  functions  Uj (x) , . . . ,un(x)  and  Vj(t),...,v  (t)  which  give  the 
Green's  function  G(x,t)  as  a finite  rank  modification  of  K(x,t)  are 
then  found  by  solving  for  the  initial  conditions,  in  terms  of  integral 
transforms  of  h(x)  to  satisfy  the  given  boundary  conditions,  it  being 
assumed  that  the  resulting  system  of  equations  has  a unique  solution.  Of 
course,  only  selfadjoint  boundary  value  problems  give  rise  to  symmetric 
Green's  functions.  In  the  next  section,  nonsynmetric  separable  kernels 
will  be  discussed. 

For  nonl inear  boundary  value  problems, 

(5.53)  h(x)  = f(x,y(x))  , 

the  use  of  the  Green's  function  G(x,t)  leads  to  a Hammerstein  integral 
equation 

1 

(5.54)  y(x)  - / G(x,t)f(t,y(t) )dt  = 0 . 

0 
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If  G(x,t)  is  of  the  form  (4.1),  then  (5.54)  is  equivalent  to  the  non- 
linear system 


(5.55) 


r x 

y(x)  - / K(x,t)f(t,y(t))dt  = au(x)  , 
0 

1 

a = / v(t)f(t,y(t))dt  . 


V 0 

The  first  equation  of  (5.55)  is  a nonlinear  Volterra  integral  equation. 
If  this  can  be  solved  for 


(5.56)  y(x)  = y(x;a)  , 

then  substitution  into  the  second  equation  gives  the  single  nonlinear 
scalar  equation 

1 

(5.57)  a = / v(t)f(t,y(t;a))dt  = <|>(a) 

0 


for  a.  In  general,  this  procedure  for  solution  cannot  be  carried  out 
explicitly,  as  in  the  linear  case,  so  various  approximation  methods, 
usually  based  on  iteration,  have  been  studied  for  this  problem.  In 
particular,  it  may  happen  that  a = y'(0),  in  which  case  iterative 
determination  of  a and  hence  y(x)  is  sometimes  referred  to  as  a 
"shooting"  method  [10].  If  G(x,t)  is  of  the  form  (4.30),  then  a similar 
construction  leads  to  a single  nonlinear  Volterra  integral  equation  for 
y(x)  = y(x;<jj,...«a  ) and  the  nonlinear  scalar  system 
1 

(5.58)  a.  = / v^(t)f(t,y(t;a1 an))dt,  1 = 1,2,. ...n  , 

0 

for  ai ,a2 <*n.  Further  discussion  of  nonlinear  boundary  value  problems 

is  outside  the  scope  of  this  paper. 
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6.  The  nonsymmetric  case.  In  certain  problems,  one  may  have  to  deal 
with  separable  kernels  which  are  not  symmetric,  such  as  finite  linear 
combinations  of  kernels  of  the  form  r(x)G(x,t)w( t) , where  G(x,t) 
is  a symmetric  simple  separable  kernel  (5.1).  As  a prototype  of  separable 
kernels  in  the  general  case,  consider  the  simple  kernel 


(6.1) 


G(x,t)  = 


u(t)v(x),  0 < t < x < 1 , 


P(x)q(t), 


subject  to  the  condition 


(6.2) 


u(x)v(x)  = p(x)q(x) , 


0 < x < t < 1 , 


0 < x < 1 , 


which  insures  that  the  traces  of  G(x,t)  and  its  associated  kernels 
are  uniquely  defined.  Formulas  will  now  be  developed  for  the  Fredholm 
determinant  a( x ) and  the  resolvent  kernel  r(x,t;X)  of  the  simple 
separable  kernel  (6.1).  The  general  separable  kernel,  which  is  a finite 
linear  combination  of  kernels  of  the  form  (6.1),  can  then  be  handled  by 
the  technique  given  in  §4,  or  the  method  to  be  discussed  in  §7. 

Writing  the  Fredholm  integral  equation  (4.2)  as  the  system 


(6.3) 


where 

(6.4) 


r y(x)  - a / «(x,t)y(t)dt  = f(x)  + Acp(x)  , 

J 0 


c = / q(t)y(t)dt  , 
0 


K(x,t)  = u(t)v(x)  - p(x)q( t)  , 


one  obtains  as  before 


(6.5) 


1 1 

a(x ) = 1 - A / q(x)P(x)dx  = 1 - A / Q(t)p(t)dt  , 
0 0 


with  the  functions  P(x),  Q(t)  given  by 


(6.6) 


r * 

P(x)  = p(x)  + A / R(x,t;A)p(t)dt  , 
0 


1 

Q(t)  = q( t ) + A / q(x)R(x,t;A)dx  , 
t 


in  terms  of  the  Vol terra  resolvent  kernel  R(x,t;A)  of  K(x,t).  If 
a(a)  i 0,  then  equation  (4.2)  has  the  unique  solution 

(6.7)  y(x)  = f(x)  + A / R(x,t;A)f(t)dt  + J1  P(x)Q(t)f(t)dt  , 

0 0 

from  which  it  follows  that  the  Fredholm  resolvent  kernel  of  the  kernel  (6 


R(x,t;>)  + ^ 


(6.8) 


r(x,t;A)  = 


P(x 


A(  A 


0 < t i x < 1 , 


0 < x < t < 1 . 


Alternatively,  one  may  write 


1 


(6.9) 


y(x)  + A / K(x,t)y(t)dt  = f(x)  + Adv(x)  , 
x 

1 

d = / u(t)y(t)dt  . 


This  leads  to  the  expressions 
(6.10) 


1 1 

A(  A ) = 1 - A / u(x)V(x)dx  « 1 - A / U(t)v(t)dt 
0 0 


for  the  Fredholm  determinant  of  G(x,t)  in  terms  of  the  functions 


(6.11) 


r i 

V(x)  * v(x)  - A / R( x , t; A ) v( t)dt  , 
x 


U(t)  = u(t)  - A / u(x)R(x,t;A)dx  , 
0 
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1)  is 


and,  if  a( a)  f 0,  to  the  unique  solution 


(6.12)  y(x)  = f ( x)  - A / R(x,t;x)f(t)dt  + / V(x)U(t)f(t)dt 


of  (4.2).  Thus, 


(6.13) 


r( x , t;  a ) = 


U(t)V(x) 

a(a)  ’ 


0 < t < x < 1 , 

0 < x < t < 1 , 


provided  that  the  Fredholm  determinant  does  not  vanish.  In  (6 . 1 1 )-(6. 1 3) , 
use  is  made  of  the  extension  of  the  function  R(x,t;x)  defined  by  (2.2) 
to  the  triangle  0 < x < t < 1. 

Comparison  of  (6.13)  and  (6.8)  yields 
(6.14)  A(A)R(x,t;x)  = U( t)V(x)  - P(x)Q(t),  0 < x,t  < 1 , 


and  finally, 

(6.15)  r(x,t;A)  = 

Theorem  6.1.  If  A(x)  f o,  then  the  simple  separable  kernel  (6.15) 
is  the  Fredholm  resolvent  kernel  of  the  simple  separable  kernel  (6.1). 

If  a(a*)  = 0,  then  it  follows  from  (6.14)  that 

(6.16)  U*(t)V*(x)  = P*(x)Q*(t)  . 


where 

(6.17) 


U*(t)  = U(t;x*) , Q*(t)  = Q( t ; A* ) , 
V*(x)  = V(x; A*),  P*(x)  = P(x;a*)  , 


are  obtained  from  (6.6)  and  (6.11)  with  A = A . It  follows  from  (6.16) 


that  V (x)  and  P (x)  are  linearly  dependent,  and  either  may  be  taken 
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as  a right  eigenfunction  of  G(x,t)  corresponding  to  the  eigenvalue 
★ ★ ★ 

x ; similarly,  U (t)  and  Q (t)  are  proportional  and  furnish  a correspond- 
ing left  eigenfunction  of  the  kernel  (6.1). 

Of  course,  the  functions  U(t),  V(x),  P(x),  Q(t)  may  be  found 
by  solving  the  appropriate  Volterra  integral  equations  if  it  is  desired 
to  avoid  the  explicit  computation  of  R(x,t;A). 


7.  An  alternative  computational  method.  Suppose  that  G(x,t)  is  a 

Fredholm  or  Volterra  kernel  with  known  resolvent  kernel  r(x,t;A),  and 

it  is  desired  to  construct  the  resolvent  kernel  r (x,t;A)  of  the  modified 

n 

kernel 


(7.1)  G ( x, t)  = G(x,t)  + l u.(x)v.(t)  . 

n j=1  J J 

In  the  special  case  that  G(x,t)  = K(x,t)  is  a Volterra  kernel,  one  form 
of  the  solution  of  this  problem  is  given  by  (4.44),  which  requires  the 
inversion  of  an  n x n matrix.  It  is  also  possible  to  obtain  the  resolvent 
kernel  by  a step-by-step  process  similar  to  an  elimination  method  for 
matrix  inversion.  First  of  all,  consider  the  rank  one  modification 

(7.2)  G^ (x,t)  = G(x,t)  + u1(x)v1(t)  . 

By  applying  the  method  of  §4  (solving  equation  (4.2)  with  the  kernel  (7.2)1, 
one  obtains 


1 1 

(7.3)  a,(a)  = 1 - a / v,(x)U,(x)dx  = 1 - A / V,(t)u,(t)dt 

i 0 0 


for  the  Fredholm  determinant  of  G^(x,t),  where 


(7.4) 


r i 

( x ) = ( x ) + A / r(x,t;A)u1(t)dt 

0 


1 


V,(t)  = v.(x)  + A / v, (x)r(x,t;A)dx  . 
i o 


Theorem  7.1.  If  a,  (a)  f 0,  then  the  resolvent  kernel  r^(x,t;A) 


of  the  rank  one  modification  Gj(x,t)  of  G(x,t)  is 


(7.5) 


U, (x)V, (t) 

r i (x,t;A)  - i ( x , t , a ) + ~ [ \ ) 
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which  is  a rank  one  modification  of  the  resolvent  kernel  r(x,t;>). 

This  result  has  been  exploited  a number  of  times  previously;  see, 
for  example,  (4.15),  (4.16),  (4.21),  (4.22),  (5.36),  (6.8),  and  (6.13). 
It  is  analogous  to  the  Sherman-Morrison-Woodbury  formula  for  finite 
matrices  [9,  pp.  123-124].  Now,  instead  of  (7.1),  one  may  consider  the 
sequence  of  kernels 


(7.6) 


Gq( x , t)  = G(x,t) 


Gk(x,t)  = G|(_1(x,t)  + uk(x)vk(t),  k = 1 ,2, . . . ,n  , 
each  of  which  is  a rank  one  modification  of  the  previous  kernel.  Taking 


(7.7) 


Tg(x,t;x)  = r(x,t;A)  , 


one  may  construct  the  corresponding  sequence  of  resolvent  kernels 


(7.8) 


rk(x,t;A)  = rkl(x,t;A)  + 


Uk(x)Vk(t) 


k = 1,2,. ..,n  , 


where 


r 


(7.9) 


and 


(7.10) 


Uk(x)  = uk(x)  + A / rk1(x,t;A)uk(t)dt  , 
1 

Vk(t)  = vk(t)  + A / vk(x)r|(_1(x,t;A)dx  , 


1 1 

Ak( A ) = 1 - A / vk(x)U.(x)dx  = 1 - A / V.  (t)u.  (t)dt  , 
0 0 


and  obtain  r (x,t;A),  provided,  of  course,  that 


(7.11) 


Ak( A ) f 0,  k = 1 ,2,. . . ,n  . 
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8.  Numerical  implications.  For  a Fredholm  integral  equation  with  a 
separable  kernel  G(x,t),  it  has  been  shown  that  methods  appropriate  to 
Volterra  integral  equations  can  be  used  to  obtain  the  Fredholm  determinant, 
resolvent  kernel,  and  eigenfunctions  of  G(x,t).  As  better  estimates  are 
available  for  convergence  of  the  resulting  expansions  than  in  the  general 
Fredholm  case,  effective  analytic  or  approximate  computations  can  be  carried 
out.  Although  explicit  formulas  can  be  obtained  only  for  very  simple 
problems,  it  may  be  that  the  Volterra  kernel  K(x,t)  is  a polynomial  or 
other  simple  function  of  x and  t,  in  which  case  a computer  can  be 
programmed  to  find  the  coefficients  in  the  expansion  (2.2)  of  the  resolvent 
kernel  K(x,t;x)  to  obtain  any  desired  degree  of  accuracy.  Eigenvalues 
of  G(x,t)  can  also  be  obtained  by  computing  zeros  of  the  entire  function 
a(a).  This  is  a more  difficult  problem,  but  can  once  again  be  done  with 
any  desired  accuracy. 

If  the  analytic  or  semi-analytic  approach  appears  fruitless  or 
uneconomical,  then  strictly  numerical  methods  based  on  numerical  integration 
may  be  used.  For  the  Volterra  integral  equations  considered,  these  lead 
to  lower  triangular  systems  of  equations  which  can  be  solved  quickly 
and  accurately,  even  if  large  [14,15].  As  numerical  integration  is  usually 
much  more  accurate  than  numerical  di fferentiation,  this  technique  offers 
an  alternative  to  finite-difference  methods  for  two-point  boundary  value 
problems  for  ordinary  differential  operators. 
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kernels  of  this  class  by  solving  related  Volterra  Integral  equations  and  finite 
linear  algebraic  systems.  Applications  to  boundary  value  problems  are  discussed 
and  explicit  formulas  are  given  for  a simple  example.  Analytic  and  numerical 


or  more  general  problems  are  indicated4 
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